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1. Introduction 

The covariant approach to bound-state calculations provides a very powerful tool for the study 
of relativistic two- and three-body bound states. When applied to hadrons, they rely upon the 
knowledge of QCD Green's functions. The intrinsic sophistication of these equations, however, 
presents several difficulties in practical applications due to the complexity of taking such calcula- 
tions beyond the simplest truncation known as Rainbow-Ladder (RL). 

It is only in recent years that a unified approach has been reached for the study of mesons 
and baryons, with quark-diquark calculations [jl|-|5|] surpassed by their more intricate three-body 
description @-||]. At the same time however, significant technical progress has been made in 
the covariant treatment of mesons beyond the RL truncation [|9|-[fl]]. In the case of baryons, these 
calculations are, due to technical difficulties, so far restricted to particles composed of quarks of the 
same mass. Therefore it is not feasible for the moment to calculate the plethora of baryon masses, 
as is done in other approaches such as constituent quark models (see, e.g. fli2| ] and references 
therein). Nevertheless, one of the main goals of the approach must be to identify the relevance 
of the different quantum-field-fheoretical interaction terms in bound-state phenomena, and in this 
respect it represents an excellent tool that is complementary to lattice QCD. 

So far, the calculation requires to model the quark-gluon interaction via an effective inter- 
action. However, it is interesting to note that studies within RL have been dominated by just one 
effective interaction, known as the Maris-Tandy model [13, 1J]. On the one hand, this dominance is 



well-earned since this ansatz performs very well indeed. On the other hand, with the swift improve- 
ment in our knowledge of QCD Green's functions from both lattice and functional approaches, it is 
possible to define different effective interactions which, presumably, capture more faithfully some 
of QCD's features. Note that this is not to say that they will perform better phenomenologically . 
Therefore, before the daunting challenge of bringing covariant baryon studies beyond RL (such as 
the inclusion of pion-cloud effects and three-body forces) or drawing conclusions about the role 
of these different contributions, one should thoroughly investigate the model-independent features 
within a given truncation scheme. 

2. Framework 

2.1 Bound-state equations 

The covariant description of bound-state equations begins with Dyson's equation for the con- 
nected and amputated «-quark scattering matrix T: 

T = K + KGqT , (2.1) 

where Go is the disconnected ra-quark dressed propagator and K is the ra-body kernel that details 
the interaction between the quarks that constitute the bound state. 

It is this interaction kernel K that proves to be challenging. For mesons we must respect 
chiral symmetry in order to realise the pion as a (pseudo)-Goldstone boson. This gives a precise 
relationship between the truncation of the quark-propagator Dyson-Schwinger equation (DSE) and 
the truncation of the 2-body kernel here. The simplest example of a symmetry-preserving kernel is 



that of Rainbow-Ladder (RL), though more sophisticated kernels have been studied M- 1 1 ] 
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2.2 Quark propagator and DSE 

The first ingredient for covariant bound-state studies is the dressed quark propagator, 

S- 1 (p)=A(p 2 )(i/ + M(p 2 )) , (2.2) 

where 1 /A(p 2 ) is the quark wave-function renormalisation and M(p 2 ) is the quark mass function. 
These scalar dressing functions are obtained as solutions to the quark DSE, 

S- 1 (p)=Z 2 S i +g 2 Z lf ( ^LfS(k)r v (k,p)D fiV (q) , (2.3) 

J (271) 

where q = k — p is the momentum of the exchanged gluon. Here Sq 1 (p) represents the bare inverse 
quark propagator, obtained from Eq. ( |2.2| ) by setting A(p 2 ) = 1 and M(p 2 ) = itiq. This bare mass is 
related to the renormalised one via Z2/W0 = Z2Z m m q , with Z2 and Z m the wave-function and quark- 
mass renormalisation constants, and Z\ f that of the quark-gluon vertex. The cardinal input in the 
quark propagator DSE ( [2j| ) are the gluon propagator D^ v (q) and the quark-gluon vertex F v (k,p). 

2.3 Rainbow-Ladder 

The quark-gluon interaction that appears in DSE for the quark reads: 

Zi f ^D^(q)r v (k,p). (2.4) 

In Landau gauge, D^ v is just the transverse projector T^ v (q) = 8^ v — q^q v /q 2 multiplied by the 
scalar gluon dressing function Z(q 2 )/q 2 . The quark-gluon vertex F v (k,p) is decomposed as twelve 
Dirac covariants, of which the minimal set in Landau gauge numbers eight. The full vertex can be 
written as the sum of its bare tree-level part plus a self-energy correction: F v (k,p) = Zi/y v + A v . 
The RL truncation requires that we replace the complicated structure of the quark-gluon vertex 



with the 7^ projection of the non-perturbative corrections. Hence Eq. (2.4) becomes 

Zxf^T^q) 7 ^- {Z lf + A(q 2 ))y v , (2.5) 

where now A(q 2 ) is the non-perturbative dressing of the y v part of the quark-gluon vertex, restricted 
to depend only on the exchanged gluon momentum. The remaining structures are as before. If one 
wishes to draw a distinction between the gluon dressing Z(q 2 ) and the quark-gluon dressings, it 
can be useful to define the function YyM{q 2 ) = Z\f + A(q 2 ). This distinction is very important 
because the gluon propagator is by now well-known from both Lattice studies and other functional 
approaches. However, if one wishes to take a purely phenomenological approach one can instead 
combine all scalar dressings into one effective running coupling, a e ff(q 2 ): 

Z u ^D^(q)r v (k,p)=Z 2 T^(q) y v . (2.6) 

In all cases, the Dirac structure remains the same, and Z\ follows from the Slavnov-Taylor identities 
to maintain multiplicative renormalisability. When it comes to discussing the possible impact of 
effects beyond RL, it is convenient to think in terms of the separate dressing functions Z and A. 
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The symmetry-preserving two-body kernel corresponding to this yields the 'ladder' part of 
RL. For simplicity we quote this in terms of CC e ff, 

^2-body = 4?rZ 2 Oefftf) ^ ^ yV ^ ?) 

q 

For the baryon, the three-body kernel ^ 3 - bod y is decomposed into a three-quark irreducible contri- 
bution K^° Ay and the sum of permuted two-body kernels K 2 '^ ody , with the subscript a indicating 
the spectator quark: 

^3-body = ^3-body + £ ^ ^2-body (2 g) 

Note that for mesons the two-body kernel describes quark-antiquark correlations, whilst in the 
baryon it pertains to quark-quark correlations. Motivated by the success of quark-diquark calcula- 
tions, we ignore Kf n body . 

2.4 Model interaction 

When one constructs a model interaction our first constraint is that of perturbation theory. This 
only fixes the large momentum behaviour, leaving the IR behaviour unspecified. This is typically 
chosen such that Dynamical Chiral Symmetry Breaking (DCSB) is realised. The chiral conden- 
sate, or equivalently the pion decay constant, determines the strength of chiral-symmetry breaking. 
These are fixed to their phenomenological/experimental values, respectively. With a symmetry 
preserving truncation protecting the chiral behaviour, such chiral properties are guaranteed and the 

1/2 

characteristic m n <=c m q behaviour is seen for 'small' m q . 

The vector meson mass is similarly determined, to a large degree, by the strength of breaking 
of chiral symmetry, albeit linearly with respect to the quark mass. That the interaction thus far con- 
structed also works well here is not surprising, since the vector mesons are IS states and similarly 
insensitive to L.S couplings present in the 7^ <g) 7^ interaction, that in RL are too attractive. Of 
course, meson bound-states with different quantum numbers are not so easily described. This is 
particularly true for the ground-state axial vectors which are significantly overbound here. 



In the present work we employ two model interactions. In the Maris-Tandy (MT) model [13 



140 the effective running coupling is given by 



a eff (q 2 ) = nr ] 1 x 2 e-^+ ^ V ■ / ... , y=q 2 /A 2 , (2.9) 



2ny m {\-e-y) 



f/A 2 



log[e 2 -l + (l+z) 2 ] 



Z=q 2 /A 2 



QCD > 



and features a Gaussian distribution in the infrared that provides dynamical chiral symmetry break- 
ing. It is characterized by an energy scale A = 0.74 GeV, fixed to give the pion decay constant, and 
a dimensionless parameter 77 . Many ground-state hadron observables have been found to be almost 
insensitive to the value of v\ around v\ = 1.8. The second part reproduces the one-loop running cou- 
pling at large, perturbative, momenta. It includes the anomalous dimension y m = 12/(1 LVc — 2A 7 /) 
of the quark propagator, and we use y m = 12/25, Aqcd = 0.234 GeV and A, = 1 GeV. Note that 
we also employ a Pauli-Villars like regulator with a mass scale of 200 GeV. The quark masses at 
jU = 19 GeV are 3.7, 85.2, 869 and 3750 MeV for the u/d, s, c, and b quarks, respectively. 
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The Alkof er-Fischer- Williams ( AFW) model [ 1 5 ] , on the other hand, is motivated by the desire 



to account for the ^(^-anomaly by the Kogut-Susskind mechanism. The effective coupling is 



constructed as the product of the gluon dressing Jlq ] and a model for the non-perturbative behaviour 
of the quark-gluon vertex [|i"7|], 

/ \ 2k / \ -K-l/2 / . \ -7o / \ -2<5o 

2\_af( x \ y \ ( ao + a uv x\ a uv x^ 



The four terms in parentheses are: the IR scaling of the gluon propagator; IR scaling of the quark- 
gluon vertex; logarithmic running of the gluon propagator; and the logarithmic running of the 
quark-gluon vertex. Additionally, the last two are constructed to interpolate between the IR and 
UV behaviour. The remaining terms are defined as follows: 

A Ay / 1 1 \ X = ^ 2//y ^ M ' 

A= TT7 + i +( y-ir auv=KY "\^-—i)> y=q2 ' K r (2 - U) 

and (Xq = S.9l5/Nc- Here, Aym = 0.71 GeV is the dynamically generated Yang-Mills scale, while 
A(/y = 0.5 GeV corresponds to the one-loop perturbative running. The IR scaling exponent is 
fc = 0.595353, and the one-loop anomalous dimensions are related via 1 + 7o = — 26o = \NcJm, 
with y m = 12/ (1 1/V C — 2Nf). We choose Nf = 5 active quark flavours at the renormalisation point 
/I = 19 GeV. The constant ^ = 0.968 is chosen such that a e ff runs appropriately in the UV. Finally, 
A/r = 0.42 GeV, As = 6.25, and Ag = 21.83 determine the IR properties of the quark-gluon vertex 
and are fitted such that the properties of %, K and p mesons are well reproduced. The quark masses 
at ji. = 19 GeV are 2.76, 55.3, 688 and 3410 MeV for the u/d, s, c, and b quarks, respectively. 

3. Results and discussion 

Using the techniques described in [^, [8]], we calculated the vector-meson, nucleon, and delta 1 
masses up to the bottom region using the AFW and MT effective interactions. The results are 
collected in Table [T], and their evolution with the pseudoscalar mass (or, equivalently, with the 
current-quark mass) is shown in Fig. |l|. 

Studying different interactions allows to quantify the model dependence within the rainbow- 
ladder truncation. We find that the two interaction models yield similar overall trends in the results. 
There is little deviation in the resulting hadron masses in the light-quark region, whereas the MT 
results tend to underestimate the AFW values for heavier quark masses. In all cases, the agree- 
ment with experimental and lattice data is < 10% and comparable to the frequently studied model 
dependence in the MT interaction that is induced by the parameter tj . 

We note that both AFW and MT interactions, which were designed in a different spirit, show 
a qualitatively similar behavior in the mid-momentum range around \q\ ~ 0.5 ... 1 GeV which is 
the relevant domain for dynamical chiral symmetry breaking. It is this region that provides the 
overall strength to ground-state hadron properties, whereas their features are less sensitive to the 
deep-infrared region. The same observation has also been made in recent studies exploring the 



impact of different model interactions in the light-meson sector [ ZjJ , 22 1. 



Assuming isospin symmetry, A and £1 baryons have the same structure, changing only the current-quark mass. 
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j pc = o- + 


MT 


AFW 


exp. 


rin (7r) 


0.140f 


0.139| 


0.138 


ris (K) 


0.496f 


0.497| 


0.496 




0.697 


0.686 




cc (n r ) 


2.979f 


2.980f 


2.980 


bb (rib) 


9.388f 


9.390| 


9.391 


J PC =\- 


MT 


AFW 


exp. 


rin (p) 


0.743 


0.710 


0.775 


ris (K*) 


0.942 


0.961 


0.892 


ss ((j>) 


1.075 


1.114 


1.020 


cc (J/y) 


3.163 


3.302 


3.097 


bb (T) 


9.466 


9.621 


9.460 





MT 


AFW 


exp. 




N 


0.94 


0.97 


0.94 




A 


1.26 


1.22 


1.23 




Q, 


1.72 


1.80 


1.67 






MT 


AFW 


lattice 


LPW 


Li ccc 


4.4 


4.9 


4.7 


4.9(0.25) 


&bbb 


13.7 


13.8 


14.4 


14.5(0.25) 



Table 1 : Computed meson and baryon masses (GeV) for both MT and AFW interactions, compared 
to experiment. Quantities fitted to their experimental values are indicated by a f . Since the heavy- 



Omega baryons have not been observed yet, we compare to lattice calculations pq , |19Q and a recent 
study from pNRQCD & 




Figure 1: Left panel: evolution of M(0) as well as p, N and A masses with the squared pion mass 
for MT and AFW models. Results are compared to lattice data; see ^ ^] for references. Right 
panel: vector-meson and triply-heavy omega masses for the two interaction models MT and AFW. 
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On the other hand, the impact of DCSB is reduced in the heavy-quark domain where hadron 
properties become increasingly sensitive to the heavy quark mass. We have included the value of 
the quark mass function M(p 2 = 0) in both panels of Fig. [|. Since in the heavy-quark limit the 
variation of M(p 2 ) at low momentum is small, its value is indicative of the heavy-quark constituent 
mass. Indeed we find that the spread between the AFW and MT results in the vector-meson and 
delta channels in the charm and bottom regions follows the same pattern as that of M(0), indicating 
that their properties are dominated by the features of the quark DSE rather than the details of the 
effective interaction and the structure of the qq kernel that enters the bound-state equations. 

In that respect it is interesting to speculate about possible effects beyond the RL truncation. 
Such contributions refer to corrections in the quark-gluon vertex as well as additional structures 
beyond the vector-vector interaction. Amongst others, they may consist of attractive pion-cloud 
corrections in the chiral regime, or repulsive corrections from self-interactions of the gluon. They 
have a two-fold effect: firstly, they allow for non vector-vector interactions and a more general 
momentum dependence; and secondly, they introduce a quark-mass dependence in the interaction 
due to the internal quark propagators that are coupled there. This last point is certainly of relevance 
when trying to describe such a wide range of mass scales provided by the u/d to b quarks. 

In rainbow-ladder, vertex corrections beyond y v are included in the modelling; otherwise Dy- 
namical Chiral Symmetry Breaking would be absent. Since this is done for the light-meson sector, 
it implies also that only light quark-dynamics are included. As a result, to determine how much 
room there is for beyond-RL effects, one must consider removing those vertex corrections that are 
implicitly included. Likely, this can only be achieved for the heavy-quark sector where one expects 
vertex corrections to be suppressed. This will be considered in a future work. 

4. Summary and conclusions 

We have studied the quark-mass dependence of several meson and baryon masses in the RL 
truncated Dyson-Schwinger approach. We investigated a broad range of current-quark masses from 
the light-quark domain up to the bottom region. To identify model-independent features, we em- 
ployed two different effective interactions for the quark-(anti)quark kernel. We find that both mod- 
els yield comparable results that agree with experimental and lattice data within < 10% throughout 
the quark-mass range, thereby demarcating the model sensitivity within a RL truncation. 

The impact of beyond-RL corrections, on the other hand, is easier to access in form-factor 
studies where they are needed for a correct description in the chiral and low-momentum region. 
In that respect it is desirable to perform a model comparison of nucleon and delta electromagnetic 
form factors in the covariant three-body framework. Work in this direction is in progress. 
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